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We give a complete structural characterisation of the map the positive branch of a one-way pattern 
implements. We start with the representation of the positive branch in terms of the phase map de- 
composition 0|, which is then further analysed to obtain the primary structure of the matrix M, 
representing the phase map decomposition in the computational basis. Using this approach we ob- 
tain some preliminary results on the connection between the columns structure of a given unitary and 
the angles of measurements in a pattern that implements it. We believe this work is a step forward 
towards a full characterisation of those unitaries with an efficient one-way model implementation. 



1 Introduction 

The one-way model of quantum computation has drawn considerable attention, mainly because it sug- 
gests different physical realisations of quantum computing (6) 13. In this model quantum states are 
transformed using single qubit measurements on an entangled state (called open graph state), which 
is prepared from an input state by performing controlled-Z operations on pairs of qubits, including the 
input system and auxiliary qubits prepared in the |+) = A= (|0) + 1 1)) state. Quantum measurements 
are probabilistic in general, and can drive the computation over 2" different branches, where n is the 
number of measurements. However, there exist sufficient conditions based on the structure of the graph 
state where the computation can be controlled by means of single qubit corrections, dependent on the 
previous measurement outcomes, so that the entire computation becomes deterministic |j6j |3] |2j |7). In 
such a deterministic computation, all the branches implement the same unitary map introduced by the 
positive branch (also known as the post-selected branch) which corresponds to the scenario in which 
every measurement collapses the qubit states to pre-selected states, typically \+ Uj ) = ^ (|0) + e' a, '|l)). 

We give a complete structural characterisation of the map the positive branch of a one-way pattern 
implements. The positive branch of a one-way pattern can be expressed in terms of a phase map decom- 
position R<t>P [4, 1], which we then further analyse to obtain the primary structure of the matrix M which 
represents R^>P in the computational basis. The columns of M can be written as: 

Me, = EiBj~$ 

where £,■ are complex scalars of norm one, parameuized by the measurement angles of the input qubits, 
Bi are signs matrices, depending on the geometry of underlying open graph state, and ~$ is a vector 
parametrized by the measurement angles of measured auxiliary qubits. The primary structure offers the 
following simple observations concerning the matrix M: 

• The first column is determined only by the geometry of the open graph state and the measurement 
angles of the auxiliary qubits. 

• All the entries of each column are sums of complex numbers of a fixed set, possibly differing in 
signs. 
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• The measurement angles of input qubits parametrize the global phase factors of columns of matrix 
M, which otherwise depend only on the geometry of the open graph state and the measurement 
angles of the auxiliary qubits. 

Moreover we can use this characterisation to easily prove the following simple lemma about uniform 
determinism. Recall that a pattern is called uniform deterministic if it is deterministic for all possible 
angles of measurements. 

Lemma 1 A pattern is uniformly deterministic if and only if it is deterministic for all possible choices of 
auxiliary measurement angles. 

We then proceed to meticulously dissect the B, matrices to reveal their structure given by the follow- 
ing decomposition: 

B, = JiAiSBNO-i, 

where ji is a sign, which depends on the adjacency of the input qubits, A,-,S,iV and £2, are diagonal sign 
matrices parameterised by the adjacencies of the set of input to the set of output qubits, the adjacency 
of output qubits, the adjacency of measured auxiliary qubits, and the adjacency of the set of input to 
the set of measured auxiliary qubits, respectively. B is a full sign matrix, parametrized by the adjacency 
between the set of output and the set of measured auxiliary qubits. The scalars and the matrices are given 
in terms of explicit functions on graphs, represented purely graph-theoretically, as adjacency matrices, 
and as lists of edges. These functions have group-theoretical properties, which we feel could further be 
utilised to help elucidate the following open problems: 

• Simulation of given unitaries directly in the one-way model, i.e. without reference to the circuit- 
based model. 

• Characterisation of graph states which implement the same map in the positive branch. 

• A refined characterisation of determinism. 

• Characterisation of the pointless measurement [5] which is a key element in defining new error 
correcting codes. 

2 Preliminaries 

In this section we briefly review the one-way model, and present the Phase Map Decomposition (U [T| 
of one-way patterns. A brief summary of linear algebra, and the notation used throughout this paper is 
given in the Appendix. 

The process of computation in the one-way model can be summarised in the following steps: 

1 . The setting up of n input qubits in an input state | y/) 

2. The addition of m — n auxiliary qubits, prepared in the state |+) = -j= (|0) + 1 1)) 

3. The pairwise entanglement of some qubits by means of the AZ interaction. This interaction is rep- 
resented by an open graph state, an ordered triplet (F,I,0), where T represents the entanglement 
graph (two qubits are entangled if and only if the corresponding vertices are adjacent), / is the set 
of input qubits/vertices and O is the set of output qubits/vertices which is a subset of the auxiliary 
qubits. 

4. The measurement of the input qubits and non-output auxiliary qubits (which we call pure auxiliary 
qubits) in the (X,Y) Bloch sphere plain, that is in the basis pairs {(| +«,-}, |— a,))}, parametrized 
by a set of measurement angles {oci, . . . , a m - n }. Here we use the following shorthand notation: 
I ice) = 75 (|0) ±e ia |l)) . The set O corresponds to the qubits which will not be measured. 
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Without loss of generality we assume input and output qubits are not overlapping. This is not a 
restriction, as additional auxiliary qubits can be added, which will correspond to the overlapping qubits, 
to which the quantum state of the overlapping qubits can be teleported. It can be easily shown these two 
scenarios are equivalent. As quantum measurements are generally probabilistic, the pattern implements 
a general completely positive map [8]. The scenario in which each measurement corresponds to the 
projection into the state state is called the positive branch, and the positive branch realises a linear 
transformation of the Hilbert space of the input qubits to the Hilbert space of the output qubits. The 
corresponding model is called projection-based quantum computing. 

We focus on the positive branch only, for this not to be a restriction, it will suffice that the graph T, 
defined by the underlying graph state, fulfils the graph-theoretical condition of having flow or generalised 
flow EUl, as then by means of local single qubit corrections, conditioned on sequential measurement 
outcomes, the entire quantum evolution of the system can be driven to be equal to the positive branch. 

We will choose the labelling of qubits so that the first n labels correspond to the input qubits, the 
following a = m — 2n correspond to the measured auxiliary qubits (which we will call pure auxiliaries), 
and the last n correspond to the output qubits. We have assumed that in a given one-way pattern all the 
input qubits are measured first (the first round of the computation). One could easily adapt the whole 
discussion of this paper to the scenario where there exist no input qubits or some of the pure auxiliary 
qubits are also measured in the first round by labelling such qubits among the first n labels. 

The measuring of a qubit in the {|±o;)} basis is equivalent to first locally rotating that same qubit 
by the local Z a unitary transformation, followed by a measurement in the {|±)} basis. For reference 
reasons, here we give the matrix representations of AZ and Z a in the computational basis: 



AZ: 
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Hence, since AZ and Z a are commuting, the projection-based computation process can be restated as 
follows: 

1. The setting up of n input qubits in the input state | y/) 

2. The addition ofm — n auxiliary qubits, set in the state |+) 

3. The application of local Z rotations to the input and m — 2n auxiliary qubits, corresponding to the 
measurement angles {ai , . . . , a m _„} 

4. The pairwise entanglement of some qubits by means of the AZ interaction. This interaction is rep- 
resented by an open graph state, an ordered triplet (r,/,0), where T represents the entanglement 
graph (two qubits are entangled if and only if the corresponding vertices are adjacent), / is the set 
of input qubits/vertices and O is an n qubit/vertex subset of the auxiliary qubits representing the 
output qubits 

5. The projection of the input qubits and m — 2n auxiliary qubits to the |+) state 

The first and second steps above comprise an embedding of a 2" dimensional Hilbert space to a 2 m 
dimensional Hilbert space which we will denote P (for preparation map), given explicitly as: 



P:\V)^\V)®\+) 



i(m—n) 
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The application of the m — n local rotations implements a map which we denote <!>! : 

m-n 

i=i 

ZiL denotes an ra-qubit unitary, which acts trivially on the composite subspaces of all qubits, except 
for the i th qubit, where it preforms the Z_«,. rotation. Note that this is an operator on a 2 m dimensional 
Hilbert space. We collect the entangling interactions, AZ, into the map <&2'. 

<J> 2 = n AZ < ; 

where the indexing goes across the set of unordered edges & of the graph state given by the graph T: 

£={{vi,v j }\{v i ,v j }QV(r)} 

The operator AZ,-j is an m-qubit unitary transformation, which acts trivially on the component subspaces 
of all qubits, except the composite subspace of qubits i and j, where it preforms the AZ transformation. 
We call the cumulative action of the latter two maps the Phase map and denote it <J> : 

<J> = <£> 2 <£> l 

The last step of the computation consists of projecting the first m — n qubits to the state |+), which we 
denote R (for restriction map:) 

/? = (+|®( m -")®/ 2 „ 

where l%n is the identity map on the 2" -dimensional Hilbert space. 

Now, the entire process of computation in the projection-based model is represented by 

R$>P 

and we call this representation the Phase map decomposition of a given unitary operator implemented 
in the one-way pattern @[T]|. Note that one can also derive directly a phase may decomposition for any 
unitary operators without any references to the one-way pattern H. 

3 Structural characterisation of the Phase map decomposition 

Let denote the standard computational orthonormal basis of a 2" dimensional complex Hilbert 

space. Every computational basis in this representation describes a sequence of 0-1 which is the binary 
representation of the integer value i — 1. Therefore i — 1 represented in binary, encodes the choice of 
states |0) or |1) in the component single qubit state spaces. For example, the state |3), in a four qubit 
setting, represents the state |0) |0) 1 1) |0) as (3 - 1) = (0010) 2 . 

Next we refine the expression R<t>P\i) to obtain the sttucture of the i th column of the matrix which 
represents R<t>P in the computational basis. 

Theorem 1 Let R<£>P be a phase map decomposition corresponding to a positive branch of a one-way 
pattern over m qubits, with n non-overlapping input and output qubits, a = m — 1n measured auxiliary 
qubits, with the set of measurement angles { — (X\ , — 0C n+a }. Then, the matrix M representing R<t>P is 
characterised with respect to columns by the following equality: 
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Me,- = (1) 

where 

• e,- is the i th vector of the canonical basis 

/"Ml \ 

• £ ; = I (^) •„ ) e r ] , w/f/i (•, •) denoting the symmetric dot product 




k=n+\ 

• B, ijS a matrix of signs of dimension 2" x2", which depends on the underlying graph state, and we 
call them the sign pattern matrices 

Proof. The proof is based on simple linear algebra manipulations so we put the details in the Ap- 
pendix. The main properties used are the diagonal form of both Z a and AZ in the computational basis. 
The complex phases arising form the Z a local rotations are collected in the ~$ vector and in the scalars 
Bj, and the diagonal of the <J> 2 entangling operation gets spread across the sign pattern matrices B ; . The 
proof itself presents this structure of the B, matrices (see Appendix) 

B i=Yj L b Ui-l)2("-")+(l-l)2"+j}\j)( l \ ( 2 ) 
j=l 1=1 

These properties will be used in the following section. □ 

For simplicity, in the expression for ^ as a numerical vector, we omit a normalising factor of 
2-( 2 V ! ) > along with the scaling factor (2~( m ~")) of the B, matrices as it has no bearing on the struc- 
ture we wish to present. 

A few direct consequences of Theorem Q] are easily checked: 

• The first column of M is parametrized by the measurement angles of pure auxiliary qubits only, as 
ei = l. 

• For every column i, the entries per row, are of the form 

Si(r,-$) (3) 

for some vector of signs r. So, entries of a column are sums of the same set of elements, varying 
possibly in sings only. 

• From [3] it is clear that every entry of every column is a sum of elements of the set of entries of the 
vector ~$ varying in signs, multiplied by the column's corresponding global phase factor £,-. 

As mentioned before we can also prove the following simple lemma about the uniform determinism. 



Lemma 2 A pattern is uniformly deterministic if and only if it is deterministic for all possible choices of 
auxiliary measurement angles. 
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Proof. Due to Theorem \T\ the measurement angles of the input qubits appear only as global phase 
factors of the columns of M, and these global factors £, are of norm one. Hence the choice of measure- 
ment angles of input qubits do not influence the norm of the columns. Also regardless of the measurement 
angles of the input qubits (as the product of two complex numbers of norm 1 is always norm 1), the ma- 
trix M is orthogonal since its columns are orthonormal. Therefore, uniform determinism can only depend 
on the measurement angles of the measured auxiliary qubits. □ 

The statement of Theorem Q] indicates a direct method for addressing problems of equalities of pat- 
terns, and of simulating a given unitary evolution of a quantum system in the one-way model. For the first 
problem, we have to evaluate and check the equalities of two expressions of the form of the right-hand 
side of [T] However, that entails knowing how to construct the sign pattern matrices from given graph 
states. This demands further analysis of the sign pattern matrices, which will be the topic of the next two 
sections. 

4 Graph-theoretical characterisation of sign pattern matrices 

In the proof of the Theorem [T] the matrices 5, were defined as representations of the expression 

H H b\(i-mi">-*)+(l-l)2*+j\\j)V\ 
7=1 1=1 

in the computational basis, and in that representation bi corresponds to the the I th diagonal entry of the 
matrix representation of <I>2 in that same basis. We now link the graph theoretical aspects of the graph 
state defining the pattern and the above expression. 

Recall that the AZ interaction is diagonal in the computational basis, hence the map <I>2 is diagonal 
in that basis as well. We introduce the sign parity function SP defined as 

SP(k) = (-lf 

as it visually simplifies the expressions. It was shown in Q that bi, the I th diagonal element of <1>2 is 
given by the following expression: 

bi = SPy x i x i i 
\(y)e* / 

where Xk was defined as the k th most significant digit (k th digit from the left, including leading zeroes) 
of / — 1 represented in binary, and S is an unordered list of edges of the entanglement graph state, 
represented by the graph T. It is easy to give a graph-theoretical representation for the expression for b\. 
It can be shown that the expression E(ij)e# x i x j > wnere Xi and xj are functions of I, counts the number 
of edges of a vertex-induced subgraph of the graph I, where the vertex set which induces the subgraph 
is determined by /. In order to explain how / determines a subset of vertices, we define a function which 
realises this vertex determination process by an integer: 

Definition 3 If V is a set of vertices, labelled with integers { 1 , . . . , m}, and k is an integer in { 1 , . . . , 2 m } 
then the selection function Sel is defined as follows: Sel(V,k) is a subset ofV such that for the vertex 
labelled with I (which we present in the subscript) V/ € V, V/ G Sel(V,k) holds if and only if the I th most 
significant digit of the m digit binary representation (including leading zeroes) ofk—lisl. 
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This function easily extends to any finite totally ordered set O, via an order-preserving bijection 
between O and {1,...,|0|}. Also, we will use the expression of the form a subset of S, selected by 
(the integer) k to mean precisely Sel(S,k). Using the introduced terminology, bi is the sign parity of the 
number of edges of the vertex-induced subgraph of T induced by a subset of the vertices of T, selected 
by I. 

Now, we can direct our attention to the expression [2] and state the following proposition about the 
graph-theoretical characterisation of the sign pattern matrices B, . 

Proposition 4 Let T be the a graph of a graph state, with vertices labelled by integers {1, . . . ,m}, such 
that the first n and last n correspond to input and output vertices (qubits) respectively, and the remaining 
a = m — In vertices correspond to the measured auxiliary {pure auxiliary) vertices {qubits) then every 
entry , is a sign parity of the number of edges of a vertex-induced subgraph ofT, and the inducing 

set of vertices V' depends on the triplet (i,p,q) as follows: 

V = Sel(I,i)USel(Aux,q)USel(0,p), 

where O denotes the subset of output vertices, Aux denotes the subset of pure auxiliary vertices, and I 
denotes the subset of input qubits. 

Proof. From the graph-theoretical characterisation of the diagonal elements of <J>2 the expression for 
b'jS and the binary representation of the index of b in |2] it is easy to see that the Proposition holds. □ 

Using the terminology of the selection function we can restate this proposition in the following 
fashion. The (p,q) entry of the sign pattern matrix is the sign parity of the number of edges of a vertex 
induced subgraph of T. This inducing subset is a union of subsets of the input, output and pure auxiliary 
vertices. The index of fi ; (the corresponding column of M) i selects a subset of the input vertices. The 
row p selects a subset of the output vertices. Finally, the column q selects a subset of the pure auxiliary 
vertices. 

Theorem Q] and Proposition @] could be potentially used in address the following problems. The 
equality of patterns and the simulation of a given unitary. In doing so, the essential expression we need 
to calculate is the expression [3] If we are interested in verifying the equality of two patterns, we need to 
calculate and compare the matrices of their phase map decompositions, given by the Theorem [T] This 
entails calculating the dot product of the rows of the matrices B, and the vectors ~$ . Similarly, if we are 
trying to simulate a given unitary, expressions [3] which will contain variables, as we go across all entries 
of all columns of the matrix M, will form the left-hand sides of a system of equations we will have to 
solve (the right-hand side being the entries of the given unitary). 

The dot product of rows of the sign pattern matrices and the vector ~$ is in general hard to evaluate, 
as both vectors have an exponential lengths in the number of pure auxiliary qubits. However, the vector 
~0 is represented as a Kronecker product of vectors of length 2, as it corresponds to a state space vector 
which can be represented as a tensor product of the minimal, two-dimensional component spaces. Such 
a representation contains the same number of 2-dimensional vectors, as there are measured auxiliary 
qubits, and so is efficient. The ability to represent the rows of the sign pattern matrices in such a compact 
form might assist in deriving techniques for solving and evaluating such expressions efficiently. 

Hence, in the following section we focus our attention to the structure of rows of the sign pattern 
matrices and present the decomposition theorem for the sign pattern matrices. 
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5 Decomposition of the sign pattern matrices 

If we turn our attention to any row p of any matrix B,, from Proposition [4] we can see that by selecting the 
index i (equivalently, a column of M) and a row p we have fixed a subset of the input qubits and a subset 
of output qubits, respectively. The p th row of B, is then generated by the sign parities of the numbers of 
edges of the vertex-induced subgraphs of T, where the inducing set is a union of the selected fixed sets of 
input and output vertices, and the entry of that row (the column of Bj) then selects the additional subset 
of the pure auxiliary vertices. 

Therefore, for fixed p and i, the corresponding row of B,, which we denote by r, is given entry-wise 
by the following expression: 

(r)k = SP (#E (ySel(lj)\JSel(0,p)\JSel(Aux,k) ) ) ( 4 ) 

where #E (T) denotes the number of edges of the graph T, and for a given graph T over the set of vertices 
V, and V' C V, Vy denotes a vertex-induced subgraph of the graph T induced by the set V' . In expression 
@]only the subsets of pure auxiliary vertices change as we go across the entries of r. 

The subgraph inducing vertex subset is expressed as a union of three subsets, two constant, and one 
variable. Let us denote A = Sel(0,p), B = Sel(I,i) and X = Sel(Aux,k). As we will be dealing with 
only one graph at a time, we will drop the graph designation and use the shorthand #E(A) instead of 
#E(Ta)- Also, with #E(Y o Z) we denote the number of edges joining vertices in Y with vertices in Z 
in the graph we are observing. It is then easy to see that 

#E(AUXUB) = #E{A) +#E(X)+ #E(B) + #£(B o A) + #E(A <^X)+#E(B^ X) (5) 

Equality [5] and the fact that the sign parity function is a homomorphism from additive monoid of integers 
to the multiplicative monoid of integers (SP(i + j) = SP(i)SP(j)) will give a basis for the decomposition 
of the sign parity matrices. Therefore, we can express the vector r with respect to the entries as follows: 

(r) k = SP (#E (B) )SP(#E(B f>A)) SP (#E (A)) SP (#E (A -H- X)) SP (#E {X )SP(#E(B hX)) (6) 

Note the dependencies of the factors of the right-hand side of [6] with respect to the explicit parameter k 
of r, parameter p which is the row selection of B, and parameter i itself which is the choice of the column 
of M Bj. 

1. SP (#E(B)) depends on i only, as it corresponds to a choice of the subset of output vertices. 

2. SP (#E(B A)) depends on both i and p, but is independent of k. 

3. SP(#E(A)) depends on p only as it corresponds to a choice of input vertices. 

4. SP (#E(A O X)) depends on p and k. 

5 . SP (#E (X ) ) depends on k only. 

6. SP (#E(B O X)) depends on i and k. 

We have represented the fixed row of a sign pattern matrix r by its entries. We will now represent r 
by using vector functions, defined on graphs, as that will allow for a simple characterisation of B, matrix 
entries. 

First, we note that, in the list of dependencies of factors which make up an entry of f, the first three 
are constants in k. The last three factors depend on k, and we shall represent them as components of 
values (which are vectors) of two different vector functions on graphs attain. 

We define a function on simple graphs whose set of vertices are equipped with a strict order. 
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Definition 5 Let F be a simple graph, where the set of vertices is equipped with a strict order. We define 
£P(F) to be a vector of signs of length 2' v given by the following components 

(&(F)) k = SP(#E(Sel(V,k))) 

for all k = 1, . . . ,2^1 

Since we will often be expressing the 2? function of some vertex-induced subgraph of a graph, it is 
convenient to adopt a shorthand notation. If the graph T, which we talk about is clear, and S is a subset 
of its set of vertices, then 3?{S) will be shorthand for 3^(Fs). Recall that T$ denotes the vertex-induced 
subgraph of the graph T, induced by the set of vertices S. 

The other useful function is defined on bipartite graphs. 

Definition 6 Let F be a bipartite graph with partitions V and W, where the set of vertices is equipped 
with a strict order. We define 3$r(V,W) to be a vector of signs of length 2^1, given by the following 
components 

(3S r (V,W)) k = SP(#E(VUSel(W,k))) 

fork = l,...,2 ]W] - 

Again, if the graph F is clear from context, we will omit the subscript F, and simply write 3S(V,W) 
instead of £$r(V, W). These two functions can be explicitly defined on different representations of graphs, 
and these representations have potentially useful properties. We give these properties after we have given 
the theorem about the decomposition of the sign pattern matrices. 

The row r can now be expressed (as its transpose, that is as a column) using 38 and 3? functions. As 
the goal is to represent a general column r (that is, any of the rows of any matrix B{), we introduce these 
parameters for row r - its row index p, and its sign pattern matrix denoted by i. Therefore, r p \ is now 
expressed as: 

r p i = YiC l p i c 2 p ■ (3§(Sel(0,p),Aux) S?{Aux) ®3§(Sel (/,/), Awe)) (7) 

where /, Aux and O denote the sets of input, pure auxiliary and output vertices and denotes the point- 
wise product. The order of the components corresponds to the order of factors in the entry-wise repre- 
sentation of r in [6] 

InEZ]?;- is a scalar, corresponding to the factor SP(#E(B)) in [6] So we can represent it using the 
function as 

Also, c p i is a constant scalar for every entry of a fixed row (hence depends on the row, and the choice of 
Bj), and corresponds to the expression SP(#E(B A)) and it can be represented using the 2$ function 

c Pti = {£S{Sel(I,i),0)) p 

Finally, c 2 p is a constant scalar for a fixed row, and does not depend on the choice of S,, and corresponds 
to the term SP (#E(A)). It can be represented using the function 

c 2 p = mo)) p . 

The three row-wise constant factors have been defined as components of vectors which depend on i 
or are constant. Then, by collecting the components across rows, and indexes i we can easily note the 
following deconstruction of the sign pattern matrices. 
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Theorem 2 Let V = IUAuxUO be the set of vertices of the graph T, tri-partitioned into input, auxiliary 
and output vertices. Let 

• Ai = diag{.^{Sel(I,i),0)) 

• S = diag (&>(())) 

• B= [@(Sel(0,l),Aux),. . . ,^(Sel(0,2^),Aux)] r 

• N = diag{& {Aux)) and 

• Q.i = diag(38(Sel(I,i),Aux)) 
then 

Bi = YiAiSBNO-i. 

Proof. The origin of j\ is straightforward and the A,- and S matrices are a direct consequence of the 
8? and 8$ function representations of the scalars c l p t and c 2 p given above. 

The B matrix contains the first factor in the brackets in the expression |7]in each row, which is constant 
in i, but variable in row p. 

Matrix N is the second factor in brackets in expression |7] spread across the diagonal of a diagonal 
matrix. That factor was constant in p and i and by presenting it as a diagonal matrix which multiplies 
B from the right, we achieve the pointwise multiplication of each row of B with that factor. Analogous 
reasoning is used for the matrix £2, with the difference that it is variable in i. □ 

Collecting the results of theorems [Hand [2] we get the following corollary: 

Corollary 1 Using the notation of theorems [7] and \2\ the matrix M can be represented with respect to 
columns as: 

MQi = EiYiAiSBNQ.j~$ . 

While Theorem [2] completely decomposes the sign pattern matrices, for an actual calculation of the 
presented decomposition it is useful to have the explicit forms of the functions SB and In the following 
section we present different representations of these functions, and present some of their properties which 
could be helpful in the application of the decomposition of Corollary Q] 



6 Explicit representations of the and SS functions 

The function 3$ has an elegant representation in terms of the adjacency matrix of the bipartite graph. If 
srf is the adjacency matrix of a bipartite graph, with partitions V and W, and all the labels of V precede 
the labels of W, then it is of the block form: 



c 
C T 



We now define the function S3 : {0, 1}" ->■ {-1, l} 2 " , such that 



&(bi,...,b„ 



i=\ 



SP{bi)_ 



(8) 
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It can be shown that if v is the modulo 2 sum of the columns of C, then 

£(V,W) = &(y). 

That is, the 3B function can be calculated directly from the adjacency matrix of the bipartite graph in ques- 
tion, by using the 38 function. Moreover, the 3$ function is a monomorphism from the group ({0, 1}" , ©) 
to the group ({ — 1, l} 2 ,0), where © and represent modulo 2 addition and pointwise multiplication, 
respectively. 

The 3$ representation of 38 and the monomorphism property are important as described below. Given 
the adjacency matrix of the underlying graph we can efficiently compute a polynomial number of entries 
of the matrix-vector multiplication Bl$ (Corollary Q]), even though the mere length of a row of B is 
exponential in the number of auxiliary qubits. It will suffice to use the representation Zj§ given on the 
right-hand side of [H and ~$ represented in the Kronecker product form, and use the following property 
of the scalar product on tensor spaces: 

(0x I ,(g)^)=n(x ! ^), 

i=\ (=1 i=l 

when Xj and F, are of equal dimensions. 

The monomorphism property also helps in the scenario where we want to calculate B£lf$ . Since 
both the rows of B and the diagonal of £2, are represented by the 3B functions, and hence by the 33 
functions, due to the monomoiphism property, the pointwise product of a row in B and the diagonal of 
£2; is again representable in the 3$ form, which easily reads out of the adjacency matrix, so this becomes 
efficiently solvable as well. 

However, there remains the problem of the matrix N, as what we really wish to calculate is BNQ^, 
which is represented by the & function. The 3? function results the sign parities of the number of edges 
of all subgraphs of a given graph as a binary vector. 

One way to explicitly represent it is by taking the positive part of the directed adjacency matrix of the 
given graph Y. That is, we direct the graph in an arbitrary fashion, and in its directed adjacency matrix 
(which carries 1 and —1 depending on the direction of the directed edges, of the now directed graph) 
replace all — l's with zeroes. If si is that matrix then it can easily be seen that 

mr))i = sp(^\i] 2 ,[i\ 2 )), 

where [i] 2 is the binary representation of i — 1 given as a vector. 

If n is the number of vertices, this representation takes n 2 binary digits on input, as they make up the 
si matrix. 

An alternative representation uses (Z\ binary digits in the form of an edge binary list, which we now 
define. Let E be an ordered set of pairs of vertices of T such that the label of the first vertex in a pair is 
strictly smaller than the label of the second, all in all Q) of them, and let E be ordered lexicographically 
according to edges; 

E = {( Vi ,Vj)\v k eV &i<j}. 

Then, for a given graph F,V = V (F) with £ we denote the binary vector of length (X) such that the i th 
entry of £ is 1 if the i th pair of vertices of E is adjacent in T and otherwise. We call this vector the 
edge binary list. It is easy to see that the edge binary list uniquely characterises a simple graph. If T is a 

graph, and £ = ( b\ , . . . , bin\ ) its edge binary list then 3? (r) can be explicitly given as 
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Wr)) ; 



where 



and 



with 



bi,...,b ( n 



Xi(i,k) 



X 2 (i,k) 



2fW 



12 



n(-i 



k=l 



mod 2, 



mod 2, 



/(*) 



The unappealing functions Xi and X2 can be explained more simply. Let 1 



be the k tb entry of the 



set Then X\ (i,k) is the <7 ?,! binary digit of binary represented i — 1, counting from the least significant 
digit. With the same notation X2(i,k) is the p th binary digit of binary represented i— 1, counting from the 
least significant digit. This representation, even though seems to be the least elegant has one significant 
properties. For & defined on edge binary lists, 

^:{0,1}(9^{-1,1} 2 " 

is a monomorphism from the group ^{0, 1}( 2 ),©^ to the group ({—1, l} 2 ",©) where © denotes point- 
wise modulo 2 addition, and pointwise multiplication. 

How to use this, or any other representation of the &P function to help efficiently evaluate or express 
N~$, or BNQ.^ in conjunction with the representation of 3§ remains an open question. 



7 Discussion 

We have presented a complete structural characterisation of the positive branch of a one-way pattern in 
terms of its matrix representation in the computational basis. This structure was shown to be intricate 
and complex yet admitting a high degree of regularity. While it remains unclear how to directly use 
this regularity to tackle problems such as direct simulation of unitaries in one-way model or full char- 
acterisation of pointless measurements and etc., however the proposed structure clearly emphasises the 
importance of the entanglement. Here, entanglement plays a crucial role in the mathematical structures 
which arise from mathematical descriptions of the process of quantum computation; If the pure auxiliary 
qubits are unconnected (unentangled), the matrix N, of the decomposition of Theorem |2j is the identity 
matrix. In that case, all the entries of the matrix realised by R<t>P can be quickly evaluated, once an 
open graph state and the measurement angles are given. If the pure auxiliary qubits are connected, this 
becomes an exponential task. We get a similar effect if we try to solve one restriction of the problem of 
simulating a given unitary. In this restricted problem an open graph state is given with the unitary, and it 
is promised that for a certain choice of angles, the positive branch will implement that unitary. For this 
promise problem it can be shown that it is easily and efficiently solvable if the pure auxiliary vertices of 



Vedran Dunjko & Elham Kashefi 



97 



the given graph are unconnected, for some families of graphs. Clearly, entanglement is again crucial. It 
is our belief that additional work on understanding the algebraic properties of the & function, that is, of 
the graph states represented as sign patterns, may yield efficient algorithms for some instances of hard 
open problems in the one-way model. Such solved instances can benefit the understanding of quantum 
computation in general. 
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8 Appendix 

8.1 Summary of notation 

Here we present a brief summary of the notation used throughout the paper. The algebra used is presented 
in the Dirac notation. 

Qubit states A qubit is represented by a two-dimensional complex Hilbert space, called the qubit's 
state space. A qubit state is a vector of unit length in the qubit's state space. The state space of an 
ensemble of qubits is represented by the tensor product of the component state spaces, and the state 
of an ensemble of qubits is a vector of unit length in the state space of the ensemble. With |0) and 
1 1 ) we denote unit orthonormal vectors in the state space of a qubit, and they constitute the standard 
computational basis of a qubit. |± a ) denotes a vectors parameterised by the real angle a (and the choice 
of + or — ) defined with respect to the computational basis vectors as 

|±„> = -^(|0>±^|1». 

When a = 0, we simply write |±). 
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Unitary transformations Z a denotes a family of phase shift unitary transformations, parametrized by 
the real angle a, represented in the computational basis with the following matrix: 



When the Z a rotation is applied to the i qubit of an ensemble of m qubits, the transformation of the 



Here, / denotes the identity operator on a single qubit state space. AZ denotes a unitary transformation 
on the state space of two qubits. In the computational basis it is given by the following matrix: 



Note that this operator cannot be represented as a tensor product of single qubit transformations. Hence, 
it can be used to create entangled states, which are multi-qubit states which cannot be represented as 
tensor products of single-qubit states. 

When the AZ transformation is applied to the component state subspace of the i th and /* qubit of an 
ensemble of m qubits, the transformation of the entire ensemble is denoted with AZ,j. The eigenvectors 
of the AZ,-j transformation are the vectors of the computational basis of the ensemble, with eigenvalue 
— 1 if both the i th and /* qubit are in the state 1 1) and eigenvalue 1 otherwise. 

Miscellaneous 

• e denotes the basis of the natural logarithm. 

• e,- denotes the i th vector of the canonical basis, i.e. a vector with entries everywhere, except a 1 
at the i th entry. 

• ® represents the tensor product. X m denotes the n-th tensor power of X, explicitly 



The tensor product of matrices (and also numerical vectors, as they are isomorphic to single row 
or column matrices) is called the Kroenecker and defined explicitly as follows: 

If A is an m-by-n matrix and B is a p-by-q matrix, then the Kronecker product A ® B is the block 




AZ = 



10 
10 
10 
-1 




matrix 



ai\B 



a\ n B 



A®B = 



a m \B 



• SP(n), for an integer n denotes the sign parity function defined as 



SP(n) = (-!)". 
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• X r denotes the transpose of the matrix (or vector) X. 

• (•, •) denotes the symmetric dot product; If X = [x\,... ,x„] T and Y = [yi, . . . ,y„] T are vectors, then 

(x,y) = £x m . 

• If r is a graph, and A C V (T) a subset of the vertices of T, Fa denotes the vertex-induced subgraph 
of the graph T induced by the set of vertices A. 

• If r is a graph #E(F) is the number of its vertices, i.e. #E(F) = |£ (r)|. If Fa is a subgraph of T, 
the graph designation can be dropped and #E(A) denotes #E(Fa)- If A and B are disjoint subsets 
of the vertices of the graph r #E(A f-> B) denotes the number of edges connecting the vertices in 
A to vertices in B in the graph Y. 

• © denotes the modulo 2 addition. If X = [x\ , . . . ,x n ] r and Y = [yi , . . . ,y„] T are vectors of integers, 
then 

X®Y = [x l ey u ...,x n ey n }\ 

• denotes the pointwise vector product; If X = [x\ , . . . ,x„] T and Y = [y\ , . . . ,y„] T are vectors, then 

XQY =[x l Qy u ...,x n Qy n }\ 

8.2 Proof of TheoremS] 

Let R<&P be the phase map decomposition [4 J of the positive branch of a one-way pattern over m qubits, 
n of which are input, n output, and a = m — 2n are pure auxiliary qubits. Also let \i) be a vector of the 
standard computational basis. Then we can directly derive the following: 

R$>P\i) = l 2 nR$>P\i) 

= Z%i\j){j\R*r\I) 

= L%i\j){j\R^®\+r m - n ) 

= L%Ai)(i\R<s>2^{\i)®\+r m - n ) 

= ((®Li<+«j) io) i'Ii ij>(ji**2 oo® (CuM ® i+r) 

For clarity reasons we temporarily omit the row-constant scalar {(®'l = i(+a k \) 10) 

= L%i \j){j\**2 (10 ® fe'VilW) ® l+r) 

= zj=i \j)(M+f {m - n) ®/2-*2 do 8> ® i+n 
= ij) ® oi) <i>2 (io ® ® \+r n ) 

We note that <I>2 = Y$=i (l\ where bi is the I th diagonal element of the diagonal matrix <I>2, 
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= Ijy)(((+r (m -" ) ®oi)l^iO(^ij (k-)^te n + ii+«,))®i+n 

2" / 2< m -") \ 
= IU) ( I hl-m»+j]( l \U\j (\i)^(^n n +l \ + a k ))^\+f n ) (9) 

= IU) f E %-i)2» + i](^7l(l0®teVil+«,))®l+>^)J 

= Ej>>( 2 E ) %-i^</i(io®(8&ii+«*>))) do) 

2" /2< m - 2 ") \ 

= E !•/')( E fe [(i-l)2(»— )+(/-l)2»+y]( Z l (®fc=n+l !+<%>)) 

= E E ^[(/-l)2('«-«)+(/-l)2»+j]U')( Z l (®fc= n +il+oft)) 

y=i /=i 

So now we summarise the entire expression (reintroducing the omitted scalar): 

(/ n \ \ 2" 2< m_2 ") / m-n \ 

0<+«J 10 E E 6 [(,--i)2(— )+(/-i)2.^]l-/)< z l !+<*> ■ 

For simplicity, in (8) we omitted a global scaling factor of 2~( 2 r !! ), brought about by the scalar 
products where they are non-zero, and in (9) the global scaling factor 2 (2), caused by 

the product {j\+) m . The overall (omitted) scaling factor is 2~5. 

The expression ((<S>k=i (+aj) I 1 ')) * s a scalar which depends on the column i, and we denote it by £,-, 
also let 

B « = Ej=lE/=l ^[(i_l)2(m-»)+(Z-l)2"+7]l^)( Z l 

be a 2" x 2( m ~ 2 ") matrix expressed in the computational basis that depends on the choice of column i 
with entries in { — 1,1}. Finally, denote the numeric representation in the computational basis of the 
vector (<8>fcTn+i l+a*)) wrtn r > which is independ of the choice of the column. It corresponds to the 
quantum state of the auxiliary qubits after the local Z a rotations, but before the entanglement procedure. 
The entire expression can then be rewritten in matrix notation as: 

Me t = SiBj~$, 

where e, is the i th vector of the canonical basis. □ 



